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Abstract 

We show that a two-dimensional photonic crystal (PC) made from a non-magnetic dielectric is 
a left-handed material in the sense defined by Veselago. Namely, it has negative values of both 
the electric permittivity e and the magnetic permeability /i in some frequency range. This follows 
from a recently proven general theorem. The negative values of e and are found by a numerical 
simulation. Using these values we demonstrate the Veselago lens, a unique optical device predicted 
by Veselago. An approximate analytical theory is proposed to calculate the values of e and fi 
from the PC band structure. It gives the results that are close to those obtained by the numerical 
simulation. The theory explains how a non-zero magnetization arises in a non-magnetic PC. 
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About 35 years ago Victor VeselagoQ] considered theoretically propagation of the elec- 
tromagnetic waves (EMW's) in a hypothetical left-handed medium (LHM), where both the 
electric permittivity e and the magnetic permeability /i are negative in some frequency 
range. Since the square of the speed of light = Cg/eyU is positive, the EMW's propagate 
but they have unusual properties, like anomalous Doppler and Cherenkov effects, negative 
light pressure, and negative refraction (NR) at the interface with a regular medium. The 
latter property enables three-dimensional (3D) imaging. The origin of all these properties 
is the fact that the energy flux, S and the wave vector, k are antiparallel in the LHM. An 
important assumption of Veselago's theory is the absence of spatial dispersion, which means 
that both e and /i are k-independent. However both of them are functions of frequency. 
A generalization of Veselago's theory for the case of spatial dispersion would not be simple 
because in this case the expression for the Poynting vector includes an extra term j^. 

To create the LHM Veselago proposed to combine two independent different subsystems 
Q], one with a negative e and another with a negative /x. As an example, he considered a 
ferromagnetic metal at frequency close to the ferromagnetic resonance and below the plasma 
frequency of the metal. Following this line, the San Diego group [3] has observed a NR in 
an artificial composite system and attributed this effect to a two-dimensional (2D) LHM. 
Their system consists of two elements: 2D metallic photonic crystal (PC) creating negative 
e below a cutoff frequency and the so-called split-ring resonators creating negative ^ due to 
anomalous dispersion near the resonance. The experiments of the San Diego group are very 
important because they have ignited the whole field, started by Veselago about 35 years 
ago but completely forgotten since then. However they caused new hot debates concerning 
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13l | predicted and observed NR in PC's using for the prediction 



the dispersion relation of EMW's and methods of geometrical optics. However, the "effective 
refractive index" they found was, in general, a function of k. We think that it originates from 
the spatial dispersion of the dielectric constant for EMW with k deep in the Brillouin zone. 
Several authors 0, Q, Q| found that near the F point this index becomes ^-independent. 
As we show below, together with the negative group velocity this indeed is a manifestation 
of the LHM. 

The main idea of this work is based upon the theorem proved in Ref. Q|. It reads 
that the group velocity = du/dh in an isotropic medium is positive in the regular 



2 



media (RM)(e > 0, fi > 0) and negative in the LHM. This means that vectors du/dh 
and k are parallel in the RM and antiparallel in the LHM. This theorem follows from the 
thermodynamics requirement that the energy of electromagnetic fields is positive. 

Now we generalize this theorem to a 2D PC illustrated in the left bottom inset of Fig. ^ 
The PC represents a square lattice of period d, comprising cylindrical holes in a dielectric 
medium with the axes of cylinders along ^-direction. We consider propagation of BMW's 
with the electric field E in this direction (s-wave). The calculated |l5| low frequency spec- 
trum of BMW's in this PC is shown in Fig. ^ for the case when the dielectric material of 
the crystal is characterized by = 12, fim = 1- The fields in the PC are represented by 
the Bloch functions with a 2D wave vector ]<.= {kx, ky). One can see from Fig. Qthat in the 
frequency range a little below the top of the band 1 at the P point, the group velocity ^uJ/^\^ 
is negative and there are no other bands in this frequency range. It is also important that 
band 1 is not degenerate at P-point. Near the P-point the wavelength is much larger than 
the lattice period of the PC and one can introduce k-independent eik and /ij^. Since 2;-axis 
is the axis of C4 symmetry, we have a uniaxial crystal and both tensors in the principal axes 
have only two components e^z, e_L = exx = ^yy and Hzz, /^x = /Wxx = IJ^yy respectively. The 
electromagnetic energy density of the s-wave can be represented in the formj2| 

U = -^&\Et (1) 

lD7ryU_|_C(j auj 

where = ezzfJ'±- In the rest of the paper we use notation ezz = e, = Since f7 > we 
get /i(vg ■ k) > 0, where 2D group velocity = 2clk/ [d^u'^n'^) / du]. Thus, negative ■ k 
means negative fi and vise versa. For the propagating waves it also means e < 0. It follows 
that this PC is a LHM in the frequency range under study. 

Note, that the p-wave (magnetic field in z direction) is characterized by e± and fizz- 
These components might also be negative in the vicinity of the P point but within a different 
frequency range. 

We find e and fi by microscopic calculations and show that they are indeed both negative 
in the frequency range under study. For this purpose we consider a PC slab infinite in 
2;-direction and embedded in a homogeneous medium with e', /i'. A 2D point source (an 
infinite dipole along z-axis) is located in the homogeneous medium to the left of the slab. If 
we now choose e' and fi' so that the cylindrical wave in the system is the same as it would 
be in a completely homogeneous medium with e' and (though inside the slab there is a 
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microscopic field), the values of e' and /i' obtained in this way are the macroscopic e and 
yU of the slab at a given frequency. The best choice for e' and yu' at the working frequency 
oj = 0.3{2iTCo/d) is obtained by a fitting procedure and the result is shown in Fig. |21 It 
yields e = —1.2, fi = —0.096. Thus, at this frequency the PC slab is a LHM and this result 
is in accordance with the above theorem. Note, that in the case of the LHM the vectors S 
and k are in the opposite directions (see animation 3]). 

One can calculate independently the refractive index of the slab ra^ = efi from the dis- 
persion relation of the band 1. At small k it can be written in a form u'^ = uf — aclk"^. Our 
calculation of spectrum gives a = 0.94, ui = 0.314(27rco/(i). Using the alternative definition 
of namely, u^n^ = c^k"^ one gets 

If. UJ 



(2) 



Note, that n — as — uoi. At the working frequency = (0.32)^, which agrees with e 
and /i found above. 

If we choose e' = lei and u' = |u| as the parameters of the homogeneous medium surround- 

n 

ing the PC, we obtain a unique device proposed by Veselago 1] that we call the Veselago 



lens. It is shown schematically in Fig. 
is shown in Fig. H] (also see animation 



1. The computer simulation of the Veselago lens 
1611). In the above simulations the PC is cut in 



the direction [10]. We have checked that the results do not change if the PC is cut in the 
direction [11]. This confirms the absence of spatial dispersion at the working frequency. The 
size of the focus is of the order of the wavelength as predicted by the diffraction theory jl7|. 

Now we explain how the starting theorem that follows from Eq.(^ and is formulated 
in terms of macroscopic electrodynamics can be interpreted microscopically. The negative 
group velocity shown at Fig. Q results from the repulsion of bands 1 and 2, 3. Therefore a 
microscopic explanation of negative /i should be found in terms of the interaction between 
these bands. We propose a quantitative explanation by writing Maxwell's equation for the 
microscopic field e^(r) in a matrix formQ, in the basis of Bloch functions at A; = 0. 
Taking into account that near the P-point the bands 1, 2, 3, 4 are far from the other bands, 
we truncate the matrix keeping only these four bands. This approximation is analogous to 
the Kane model in semiconductor physics|20]. By diagonalizing the truncated matrix we 
find the spectra of all four bands and the correct combinations of basis functions for each 
band at finite k. We have checked that the dispersion law for all bands 1-4 as obtained from 
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the truncated model does not differ by more than 4% from the result presented in Fig. at 
kd/27T < 0.1. At small k the microscopic electric field for the band 1 has the form 

/ N ,k r f 'ik^kxck ik^kyC^ ikokik'^ck \ , , 

where r is the 2D- vector, mi, M2, Ms and are the Bloch functions at A; = 0, the positive 
parameters fco, fci have the dimensionality of a wave vector and they are integrals from these 
functions. The symmetry of the function Ui is — . The function has + Y"^ 
symmetry. The functions U2 and have the symmetries X and Y respectively. The parity 
of X and Y functions within the elementary cell is shown in the central inset of Fig. ^ 

If the material of the PC were metallic, the electric field in the bands 2 and 3 would 
create two currents in opposite directions. These currents would induce a magnetic moment 
in X and y directions. In case of the dielectric PC these currents are imaginary displacement 
currents and they do not create any magnetic moment in the modes 2 and 3. However, as 
follows from Eq. (jH)) the admixture of bands 2 and 3 to band 1 is also imaginary. Therefore 
the band 1 has a real magnetic moment. 

The magnetization can be calculated using the expression 



where P = ze2(r)[e(r) — l]/47r, z is a unit vector in z direction, e(r) = 1 inside the holes and 
e^n otherwise, the integration is performed over the unit cell, and e^(r) is given by Eq.Q. 
Only the second and the third terms from Eq.Q contribute to the magnetization. The 
result can be expressed through the integrals from U2 and M3. These functions have been 
taken from the computation of the spectrum shown at Fig. ^ The macroscopic vectors B 
and E are connected by the relation B = cok x E/cu, where E = z < ez(r) > and < ... > 
means averaging over the unit cell. Only the function M4 gives a non-zero contribution after 
the averaging. It gives the main contribution to E. Other basis functions give non-zero 
contribution only when multiplied by r-dependent terms in the expansion of the exp (zk ■ r) 
within the unit cell. Their contributions to E are of the same order in k but they are much 
smaller than the contribution from since [u^ — uji)/uji ^ 1, where uJi and uj^ are the 
frequencies of the bands 1 and 4 at A; = 0. In fact, the latter strong inequality permits 
truncation of the matrix in our analytical model. 

Now we can find /i from the definition = 1 — 47r|M|/|B|. The result is 
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Comparing Eq.Q and Eq.(jHI) one can see that the root of fi at u = ui is responsible for the 
root of n^. This means that e is non-zero and negative ai uj = ui. On the other hand, n 
changes its sign aiu = uji. Thus, the gap at > cui is due to /ie < 0. At u = 0.3(27rco/(i) one 
gets fi = —0.084 and that is very close to the value obtained above by the fitting procedure. 
Thus, we have found analytically the value of fi. 

The 3D LHM can be obtained in the same way, for example, using the PC with the cubic 
symmetry proposed by Luo et a/.^^. This PC also has a negative group velocity near the 
r point with no other bands at this frequency. 

Thus, we have shown that the LHM can be readily obtained using only one dielectric 
PC that creates both negative e and fi. We have presented the physical explanation of this 
phenomenon in a 2D case. The working frequency can be easily changed in a wide range 
by changing the lattice constant and of the PC. It is important for applications that the 
losses in a dielectric PC are much smaller than in a metallic one. We have found not only 
n but also e and /i that are necessary to describe the dynamical properties of the interfaces. 
This is especially important for the creation of the Veselago lens, which requires no reflection 
at interfaces, thus imposing two matching conditions on e and fi. By matching n only one 
gets a modification of the Veselago lens with multiple foci (Fig. Eb) 2M- realize, that 
we have not presented a method of matching the LHM with the surrounding medium. This 
matching is very easy in computational approach but it must be handled creatively on the 
way toward real devices. 

We are grateful to Serge Luryi and Valy Vardeny for helpful discussions. The work has 
been funded by the NSF grant DMR-0102964. 
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FIG. 1: The lowest bands of the 2D photonic crystal that is shown in the left bottom inset. The 
cylindrical hole radius is 0.35d. Dashed line indicates the working frequency. Central insets show 
signs of electric field in a unit cell for bands 2, 3 at /c = 0. Points M and X are at the boundary 
of the Brillouin zone in directions [1, 1] and [1,0] respectively. 
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FIG. 2: Finding the values of e and /i of PC using the propagation of EMW's from the 2D 
point source through the slab of the PC. Parameters of the homogeneous medium are e' = —1.2, 
/i' = —0.096. Black curves show the circles centered at the point source. One can see that EMW's 
in the homogeneous medium almost ignore the slab. We conclude from this that e = e' , fi = fi' . 
Vectors S and k are opposite to each other as it should be in the LHM. The computations are 
performed by the finite-element method with perfectly matched layers at the boundaries. 
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FIG. 3: 3D imaging by the lenses based upon LHM. The objects are marked by bulbs. The 
negative refraction at the interfaces is obtained from electrodynamics conditions of continuity of 
the tangential component of k and the normal component of S, and from the main LHM property 
that the direction of S is opposite to the direction of k (shown by arrows), a) The lens proposed 
by Veselago (Veselago lensjlj). In this case e' = |e|, = \fi\, the reflected waves are absent and the 
lens has only one external focus, b) Multifocal lens[^| with e/x = e'/i'. The multiple foci appear 
due to reflected waves. 
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FIG. 4: Computer simulation of the 2D Veselago lens. In this case e' = |e|, fi' = The curves 
show the circles centered in the point source to the left of the PC slab and at the focus to the right 
of the slab. Details of the simulation are similar to Fig. 2. 



11 



